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Coupling of Brownian motions and Perelman's /^-functional 

Kazumasa Kuwada*^ and Robert Philipowski* 



Abstract 

We show that on a manifold whose Riemannian metric evolves under backwards Ricci flow 
two Brownian motions can be coupled in such a way that the expectation of their normalized 
^-distance is non-increasing. As an immediate corollary we obtain a new proof of a recent 
result of Topping (J. reine angew. Math. 636 (2009), 93-122), namely that the normalized 
/^-transportation cost between two solutions of the heat equation is non-increasing as well. 
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qq ; 1 Introduction 



Let M be a <i-dimensional differentiable manifold, < t\ < t 2 < T and {g(T)) Te [f lt T] a complete 
backwards Ricci flow on M, i.e. a smooth family of Riemannian metrics satisfying 



13: % = 2mc ^ (1) 



and such that (M, g(r)) is complete for all r £ [n, T]. In this situation Perelman [Tgl Section 7.1] 
(see also O Definition 7.5]) defined the /2-functional of a smooth curve 7 : [17,72] —> M (where 
f 1 < ri < r 2 < T) by 

O" 

§: £&)'■=! V^ [\i(r)\ 2 g{T) + R 9[T) (l(r)) 



n 



where R g r T ^(x) is the scalar curvature at x with respect to the metric g(r). 
q ■ Denoting by L(x, n; y, 72) the infimum of £(7) over smooth curves 7 : [71,72] — >• M satisfying 

7(ri) = x and 7(712) = y, and by 



W £ (/i,Ti;i/,T 2 ) :=inf / L(x,Ti;y,T 2 )n(dx,dy) 

n JMxM 

(the infimum is over all probability measures tt on M x M whose marginals are \i and v) the 
associated transportation cost between two probability measures \i and v on M, Topping [23] 
(see also Lott [13]) obtained the following result: 
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Theorem 1 (Theorem 1.1 in [23]). Assume that M is compact and that f\ > 0. Let u : 
[ti,T] X M — >■ R+ and v : [t2, T] x M — > R+ 6e too non-negative unit-mass solutions of the heat 
equation 

— = A g{T) u - Ru, 

where the term Ru comes from the change in time of the volume element. Then the normalized 
£ -transportation cost 

Q(t) := 2{y/fii- ^\t)W C (u(f 1 t,-)vol g ( flt ),fit-,v(f 2 t,-)volg( f2t ),f 2 t) - 2d (y/fit - \f¥J) 

between the two solutions evaluated at times f±t resp. f 2 t is a non-increasing function of t £ 

[i,r/7i]. 

By <7(r)-Browman motion, we mean the time-inhomogeneous diffusion process whose genera- 
tor is A S ( T ). As in the time-homogeneous case, the heat distribution u(t, •) vol g ( r ) is expressed as 
the law of a g(r)-Brownian motion at time r. In view of this strong relation between heat equa- 
tion and Brownian motion, it is natural to ask whether one can couple two Brownian motions 
on M in such a way that a pathwise analogue of this result involving the function 

v 2 



Q(t,x,y) : = 2{ v r r 2 l- Vnt) L(x,T 1 t;y,T 2 t) -2d{y r r 2 l- Vnt)' 

holds. The main result of this paper answers it affirmatively as follows: 

Theorem 2. Assume that M has bounded curvature tensor, i.e. 

sup | Rm 9(T ) \ g ( T )(x) < co. (2) 

xeM,re[fi,T] 

Then given any points x,y G M and any s G [1, T/ffc], there exist two coupled g{r)-Brownian 
motions {Xr)re[fis,T\ an d (XT) T e[f 2 s,T] with initial values X flS = x and Y f2S = y such that the 
process (0(i, ^fit)^f 2 t))iG[s,T/f 2 ] * s a supermartingale. In particular E[Q(t,Xf lt ,Y f2t )] is non- 
increasing. In addition, we can take them so that the map (x, y) \- > (X, Y) is measurable. 

Remark 1. Obviously, ([2]) is satisfied if M is compact. Thus it includes the case of Theorem [TJ 
In addition, there are plenty of examples of backwards Ricci flow satisfying ([2]) even when M 
is non-compact. Indeed, given a metric go on M with bounded curvature tensor, there exists a 
unique solution to the Ricci flow dtg(t) = — 2~Ric g n-\ with initial condition go satisfying ([2]) for 
a short time (see [21] for existence and [1] for uniqueness). Then the corresponding backwards 
Ricci flow is obtained by time-reversal. 

Remark 2. As shown in [12], under backwards Ricci flow g(r)-Brownian motion cannot explode. 
Hence @{t,Xt,Yt) is well-defined for all t £ [s,T/f 2 ]. This fact also ensures that u(r,-)yol g r T \ 
has unit mass whenever it does at the initial time. 

Using Theorem [2] we can prove Topping's result even in the non-compact case. 

Theorem 3. Assume that ([2]) holds. Then the same assertion as in Theorem [I] holds true for 
nonegative unit mass solutions u and v to the heat equation and the associated functional Q(t). 

Proof of Theorem using Theorem [H Fix 1 < s < t < T/f 2 , and let n be an optimal coupling of 
u(f\s, •) vol g (f ls -\ and v(f 2 s, •) vol g (f 2S y (Existence of an optimal coupling follows from J2U The- 

orem 4.1], using the obvious lower bound L(x, Ti;y, t 2 ) > g(r 2 — t x ) vai x& M,re[fi,T\ Rg(r)i x )-) 
For each (x, y) G M x M, we take coupled Brownian motions {X* ) T ^[f ls ,T] an d (Xr) T e[f2s,T] with 
initial values X^ lS = x and Y^ 2S = y as in Theorem[2j Since (x, y) h-> (X x , Y y ) is measurable, we 
can construct a coupling of two Brownian motions (X, Y) with initial distribution ir by following 
a usual manner. Then the joint distribution of X fl t and Y f2 t is a coupling of u(f\t, ■) vol g r flt \ 
and v{f 2 t,-)vol g(f2t) , so that 9(t) <E[&(t,X flt ,Yf 2 t)] < E [0(s,X flS ,Y f2S )] = &{s). D 



2 Remarks concerning related work 

The Ricci flow was introduced by Hamilton [8j . There he effectively used it to solve the Poincare 
conjecture for 3-manifolds with positive Ricci curvature. By following his approach, Perel- 
man [TH [19j [20] finally solved the Poincare conjecture (see also O [9j ED- There he used 
^-functional as a crucial tool. At the same stage, he also studied the heat equation in [TS] in 
relation with the geometry of Ricci flows. It suggests that analysing the heat equation is still an 
efficient way to investigate geometry of the underlying space even in the time-dependent metric 
case. This general principle has been confirmed in recent developments in this direction. In 
connection with the theory of optimal transportation, McCann and Topping |16j showed con- 
traction in the L 2 -Wasserstein distance for the heat equation under backwards Ricci flow on a 
compact manifold. Topping's result [23] can be regarded as an extension of it to contraction 
in the normalized /^-transportation cost (see p3] also). By taking fi — > f±, he recovered the 
monotonicity of Perelman's W-entropy, which is one of fundamental ingredients in Perelman's 
work. 

A probabilistic approach to these problems is initiated by Arnaudon, Coulibaly and Thal- 
maier. In J2J Section 4], they sharpened McCann and Topping's result [IB] to a pathwise con- 
traction in the following sense: There is a coupling (Xt,Yf)t>o of two Brownian motions starting 
from x, y E X respectively such that the <7(i)-distance between Xt and Yj is non-increasing in 
t almost surely. In their approach, probabilistic techniques based on analysis of sample paths 
made it possible to establish such a pathwise estimate. It should be mentioned that, as another 
advantage of their approach, their argument works even on non-compact M (cf. [12]). Our ap- 
proach is the same as theirs in spirit. In fact, such advantages are also inherited to our results. 
Unfortunately, we cannot expect a pathwise contraction as theirs since our problem differs in 
nature from what is studied in [1] (see Remark [7]). However, it should be noted that this new 
fact is revealed as a result of our pathwise arguments. Furthermore, we can expect that our 
approach makes it possible to employ several techniques in stochastic analysis to obtain more 
detailed behavior of 0(£,Xf lt , Yf 2i ), especially in the limit T2 — > f\, in a future development. 
Note that, from technical point of view, our method relies on the result in [IT] and it is different 
from Arnaudon, Coulibaly and Thalmaier's one. 

3 Coupling of Brownian motions in the absence of £-cut locus 

Since the proof of Theorem [2] involves some technical arguments, first we study the problem in 
the case that the ^-distance L has no singularity. More precisely, 

Assumption 1. The £-cut locus is empty. 

See subsection 15.11 or [51 [231 I2S] for the definition of £-cut locus. Under Assumption [JJ the 
following holds: 

1. For all x, y E M and all f i < T\ < T2 < T there is a unique minimizer 7™ T2 of L(x, n; y, T2) 
(existence of 7™ T2 is proved in [5j Lemma 7.27], while uniqueness follows immediately from 
the characterization of £-cut locus, see subsection 15 .1|) . 

2. The function L is globally smooth. 

Thus, in this case, we can freely use stochastic analysis on the frame bundle without taking 
any care on regularity of L. In section we present the complete proof of Theorem [2] using a 
random walk approximation (see Remark [8] for further details on the choice of our approach) . 



3.1 Construction of the coupling 

A (/(r)-Brownian motion X on M (scaled in time by the factor f\) starting at a point x G M 
at time s G [l,T/f2] can be constructed in the following way [TJ [12]: Let it : F(M) — > M 
be the frame bundle and (ej)f =1 the standard basis of R d . For each r G [^l,^] let (Hi(r))f =1 
be the associated g(7")-horizontal vector fields on J-{M) (i.e. Hi(r,u) is the g(r)-horizontal lift 
of uei). Moreover let {V a,/3 )a r=i De the canonical vertical vector fields, i.e. (V Q, ^/)(n) := 

(/(u(m))) (m = {m a p)a a =1 G GLd(M)), and let (W t )t>o be a standard Revalued 



dm aP m=w 

Brownian motion. By 9 ^ T '{M), we denote the g , (r)-orthonormal frame bundle. 

We first define a horizontal Brownian motion on J-{M) as the solution U = {Ut)t&[s,T/fi} °f 
the Stratonovich SDE 

dU t = ^2f- 1 Y,H i {T 1 t,U t )odWl-T 1 Y, ^(fit)(U t e a ,U t ep)V a P(U t )dt (3) 

«=1 «,/3=l 

with initial value U s = u E Of lS (M), and then define a scaled Brownian motion X on M as 

X t := nU t . 

Note that Xt does not move when f\ = 0. The last term in ([3]) ensures that Ut G 9 ^ Tlt '{M) 
for all £ G [s,T/fi] (see [U Proposition 1.1], [3 Proposition 1.2]), so that by Ito's formula for all 
smooth / : [s, T/fi] x M -)• K 

d/(i,l t ) = ^(t,l t )dt + v / 2^^(^e i )/(t,X t )dW t i + fiA g(Tlt) /(t ) X t )di. 

Let us define (^t)tg[tis,ti by -^"rit := Xt- Then X T becomes a g(r)-Brownian motion when 
fi >0. 

Remark 3. Intuitively, it might be helpful to think that X T lives in (M,g(r)), or Xt lives in 
(M, g(fit)). The same is true for Y and Y which will be defined below. Similarly, for all curves 
7 : [ r ii T2] — > M in this paper, we can naturally regard 7(7") as in (M, g(r)). 

We now want to construct a second scaled Brownian motion Y on A/ in such a way that 
its infinitesimal increments dYt are "space-time parallel" to those of X along the minimal C- 
geodesic (namely, the minimizer of L) from (Xt,f\t) to (Yt,f2t). To make this idea precise, we 
first define the notion of space-time parallel vector field: 

Definition 1 (space-time parallel vector field). Let f\ < T\ < T2 < T and 7 : [ti,T2] — > M be a 
smooth curve. We say that a vector field Z along 7 is space-time parallel if 

V^ ) ) Z(r) = -Ricf (r) (Z(r)) (4) 

holds for all r G [ti,T2]. Here \7 9 ^ T ' stands for the covariant derivative associated with the 
<7(r)-Levi-Civita connection and Ric , , is defined by regarding the g(r)-Ricci curvature as a 
(l,l)-tensor. 

Remark 4. Since (j4|) is a linear first-order ODE, for any £ G T^r Tl \M there exists a unique 
space-time parallel vector field Z along 7 with Z{j\) = £. 



Remark 5. Whenever Z and Z' are space-time parallel vector fields along a curve 7, their 
<7(r)-inner product is constant in r: 

±(Z(t),Z'(t)) 9{t) = ^(r)(Z(r),Z'(r)) + (V^ r ]z(r), Z'(r)) g{r) + (Z(r), V^]z'(r)) g{r) 

= 2Ric giT) (Z(T),Z'(T))-m CgiT) (Z(T),Z'(T))-mc g[T) (Z(r),Z'(T)) 
= 0. 

Remark 6. The emergence of the Ricci curvature in Q is based on the Ricci flow equation 
(PQ). Indeed, we can generalize the notion of space-time parallel transport even in the absence 
of (HJ with keeping the property in the last remark. This would be a natural extension in the 
sense that it coincides with the usual parallel transport when g(r) is constant in r. On the other 
hand, it is convenient to define it as (JU) for later use in this paper. 

Definition 2 (space-time parallel transport). For x,y £ M and f\ < t\ < t<i < T, we define a 
map rri^y 2 : T X M — > T y M as follows: mQy 2 (£) := Zfo), where Z is the unique space-time paral- 
lel vector field along 7^ T2 with Z(t\) = £. By Remark[5j m x 1 y T2 is an isometry from (T X M, g{T\)) 
to (T y M,g(T2)). In addition, it smoothly depends on x,T\,y,T2 under Assumption [TJ 

We now define a second horizontal scaled Brownian motion V = (Vt)tg[ s T/n] on J"(M) as 
the solution of 

d d rj 

dV t = V2f^Y, H i^ f ^ Vu^t) o dWi -f 2 J2 ^(T2t)(Vte a , V t ep)V a P(V t )dt 

i=l a,j3=l 

with initial value V s = v G Oy (M), and we set Y t := nVt- Here H*(u,ti;v,T2) is the 
<7(T2)-horizontal lift of ve*(u,Ti; v,T2), where 

e*(u,n;v,T 2 ) := v^rrQ^uei. 

As we did for X, let us define (Y T ) T€ r f2S ^ T i by Y f2 t := Yj to make Y a g(r)-Brownian motion. 
From theoretical point of view, it seems to be natural to work with (X T ,Y T ) (see Remark [3j). 
However, for technical simplicity, we will prefer to work with (Xt,Yt) instead in the sequel. 

3.2 Proof of Theorem [2] in the absence of £-cut locus 

Our argument in this section is based on the following Ito's formula for (JQ, Yj). 
Lemma 1. Let f be a smooth function on [s,T/f2J x M x M. Then 

df(t, X t ,Y t ) = ^(t, X t ,Y t )dt + £ [y/vRUtet © V^Vte*] f(t, X u Y t )dW\ 

i=l 
d 

+ ^ Hess 9(nt)©9(f 2 t) f\( t x t Y t ) W^U t ei © Vf2~V t e*, VnU t ei © V^Vte*) dt. 



s s(nt)es(f 2 i)/|( t)j sf t ,Y t ) 

Here the Hessian of f is taken with respect to the product metric g(f\t) © g(f2t), e\ stands for 
e*(Ut,fit;Vt,f2t), and for tangent vectors £1 € T X M , £2 £ T y M we write £1 ©£2 := (£1,^2) £ 
T (X)V) {MxM). 



Proof. Ito's formula applied to a smooth function / on [s, T/t 2 ] X J-(M) x J-(M) gives 
df(t,U t ,V t ) 

= ^(t, U t , V t )dt + £J [v^fliCfit, U t ) © V^miUunt; Vt,f 2 t)] f(t, U t , V t )dWi 



i=l 



2 _ 



d 



a,B=l 



fi-^(nt)(U t e a ,U t e )V^(U t )®f 2 -^(f 2 t)(V t e a ,V t e )V^(V t ) 



f(t,U t ,Vt)dt. 



The claim follows by choosing f(t,u,v) := f(t,iru,irv) because the function considered here is 
constant in the vertical direction so that the term involving V a P f vanishes. □ 

Let A(t,x,y) := L(x,fit;y,f2t). In order to apply Lemma Q] to the function G we need the 
following proposition whose proof is given in the next section: 

Proposition 1. Take x,y £ M, u € Ol (M) and v € Oy (M). Let 7 be a minimizer 
of L(x,fit;y,f 2 t). Assume that (x , fit; y , f 2 t) is not in the C-cut locus. Set £3 := \fr\uei © 
\fr 2 ve\(u^ fit;v,f 2 t). Then 

f)\ 1 f f2t / "3 

^(t,x,y) = -^ r 3 /^-^ (T) ( 7 (r)) - A g(r)J R ff(r) ( 7 (r)) - 2|Ric ff(r) g (r) ( 7 (r)) 

" ^rl7(r)|J (T) +2Ric ff(r) ( 7 (T),7(r)))dr, (5) 



YJ Hess 



3(7-1)09(^-2) 



A 



i=X 



(t,x,y) 



(&,& 



< 



djr~ 



T=T2t 



1 f T2t 



+ - I r 3 / 2 (2|Ric 9(T) | 3(r) (7(r)) + A 9(r)J R 5(r) ( 7 (T)) 

T=T\t J T\t 

- -R g[T) ( 7 (r)) - 2Ric fl(T) (7(r),7(r)))dT (6) 



and consequently 
dA 



01 



(t, x, y)+ Y^ H ess 9(ri)e!? ( T2 ) A| x (&, £ 



t=l 



< 



r 2 - v 7 ^) - Y t [ 2 v 7 ^ (i? ff ( T )(7(^)) + |7( 



v 7 * 

d / /^- / ^ N 1 . 
-y| (VT2 - Vn) - —A{t,x,y). 



T)\~ g{T) ] dr 



2t 



The proof of Theorem is now achieved under Assumption [T] by combining Lemma Q] and 
Proposition [1} 

Proof of Theorem^ under AssumptionUl Lemma[6]below ensures that 0(i, Xt, Yt) is integrable. 
Thus it suffices to show that the bounded variation part of @(t,Xt,Yt) is nonpositive. By 



Lemma [H 
iB(t,X,,Y t ) 



£(.,*.«) 



J 



+ Y, Hess fi 



e 



s 9(nt)eff(r 2 t) KJ \(t,Xt,Y t ) 

i=l 
d 

+ J] [V2fiU t ei © v^^tej 
i=i 

For the bounded variation part we obtain 



TiUtet © yMVte*, VnU t e,i © vWte- 

e(t,x t ,y t )dw?. 



r// 



f(t,x t ,y) 



^A(t,x t ,y) + 2 (v^t- Vn*) ^(t,l t ,y) - 2d(V^- v 7 ^") 2 



v/7 



( 9/ 



and 



y Hess 



*=l 



9(rit)©9(f2t) |(t,X t ,y t ) 



^t^ © VnV t e*, VnUta © V^eJ 



2 (v^ - V^<) X] Hess 



A 



i=l 



■s^i^es^O^I^Xt.Yi) 



/ n&te < © V^Vte^ , VnU t ei © ^V t e* ) . 



Thus, by Proposition HJ 

dO d 

— (i,X t ,y 4 ) + ^Hess 



2=1 

<2(y/^t-y/f^t) 
= 0. 



s(TLt)ffis(^t) |(t,x t ,y t ) 



-J=(^-Vn)-^A(t,x t ,y t ) 



+ Vn y n A(t, X t , Y t ) -2d{vr 2 - Vf~i) 2 



Vt 



Hence @{t,Xt,Yt) is indeed a supermartingale. 



^nUtei © Vr^V t e*, VnU t ei © v^ej 



D 



Remark 7. Unlike the case in pQ, the pathwise contraction of Q(t,Xt,Yt) is no longer true 
in our case. In other words, the martingale part of Q(t,Xt,Yf) does not vanish. We will see 
it in the following. The minimal £-geodesic 7 = 72 / 7 " 2 OI " L(x,Ti;y,T2) satisfies the £-geodesic 
equation 

V^7(t) = \V 9{T) R 9{T) - 2Ric# r) (7(r)) - ±j(r) (7) 

(see O Corollary 7.19]). Thus the first variation formula (see Lemma 7.15]) yields 



2f 1 U t e i © V2f2V t e*A{t,X t ,Y t ) = V2tf 2 (V t e*,j(f 2 t)) g{f2t) - V2tf 1 (U t e i ,j(f 1 t)) g{ 



(fit)- 



(8) 



One obstruction to pathwise contraction is on the difference of time-scalings t\ and t 2 . In 
addition, by (J7J), sfr^ij) is not space-time parallel to 7 in general (cf. Remark [5]). 



4 Proof of Proposition [T] 

In this section, we write t\ := fit and r 2 := f 2 £. We assume r 2 < T. For simplicity of notations, 
we abbreviate the dependency on the metric g(r) of several geometric quantities such as Ric, R, 
the inner product (•,•), the covariant derivative V etc. when our choice of r is obvious. For this 
abbreviation, we will think that 7(7") is in (M,g(r)) and j(r) is in (T~r T \M, g(r)). Note that, 
when 77 = 0, lim,-^ -v/r7( T ) exists while lim T j,o |7(t)| = 00. In any case, i/^'H 7 ")! i s bounded 
(see dMD). 

We first compute the time derivative of A. When 77 > 0, by |23(. Formulas (A. 4) and (A. 5)] 
we have 



dL 



(x,T 1 ;y,T 2 ) = -y/n(R g{Tl) (x) - |7(ti)| 2 ) , 



dr\ 

dL 

— (x,T 1 ;y,T 2 ) = Vn {R g ( T2 )(y) ~ \i( T 2)\ 2 ) 



so that 



<9A . _ dL _ dL 

-Q7{t, x, y) = T ig^r{x, n;y, r 2 ) + t 2 — (re, n;y, r 2 ) 



1 



.3/2 



.3/2 



- [r^ (i?( 7 (r 2 )) - | 7 (r 2 )| 2 ) - r{' A (Rfrin)) ~ \i(n)[ 



(9) 



Thus the integration-by-parts yields, 
d\ 3 f T2 



+ 



\ r T 2 



3/2 k(TW) + V 7 (r)%(r)) 



n 



V<9r 



2(V^ (t) 7(t),7(t)) - 2Ric( 7 (r), 7 (r)) dr. 



Note that we have 



^ = -A J R-2|Ric| 2 
or 



(10) 



(11) 



(see e.g. [221 Proposition 2.5.4]). Since 7 satisfies the £-geodesic equation (J7J), by substituting 
(J7|) and (|lip into ()10p . we obtain (|5|). Note that the derivation of Q and (|10p is still valid even 
when 77 = because of the remark at the beginning of this section. Thus ([5]) holds when 77 = 0, 
too. 

In order to estimate ^«=i ^- ess g(ri)®g(T 2 ) ^\u \ (£ij£i) we begin with the second variation 
formula for the /^-functional: 

Lemma 2 (Second variation formula; [5j Lemma 7.37]). Let V : (— e,e) x [77, r 2 ] — > M be a 
variation 0/7, S(s,t) := 9 s r(s,r), and Z(r) := 3 s r(0, r) £/ie variation field ofT. Then 



d?_ 

ds 2 



£(r s )=2v9(7(r),V z(T) 5(0,r)) 



s=0 



'^)| 2 
T 



r=T2 

T=Tl 

T2 



2V7 7 Ric(Z(r),Z(r))|; 



T=T2 

n 



>/rfl"(7(r),Z(r))dT 



2^/7 



V^ (T) Z(r) + Ric#(Z(r))-^Z(r; 



dr, 



(12) 



where 

H(j(r), Z(t)) := -2^(Z(r), Z(r)) - Hess R(Z(r), Z(r)) + 2 |Ric#(Z(r))| 2 

- - Ric(Z(r), Z(t)) - 2 Rm(Z(r),7(r),7(r), Z(r)) 
r 

- 4(V^ (T) Ric)(Z(r), Z{t)) + 4(V Z(T) Ric)(7(r), Z(r)). (13) 

In [5] this lemma is only proved in the case t\ = and Z{j\) = 0. However, the proof given 
there can be easily adapted to the slightly more general case needed here. 

Corollary 1 (see Lemma 7.39] for a similar statement). If the variation field Z is of the 
form 

Z{r) = ^Z*(r) (14) 

with a space-time parallel field Z* satisfying |Z*(r)| = 1, then 
d 2 



ds 2 



C(T S ) = 2^ < 7 (r), V z(T) 5(0,r)) ff(T) |^ - 2 ^Ric(Z(r),Z(r))|^ 



- r ^H(j(r),Z(T))dr + 
Proof. Since Z* is space-time parallel, Z satisfies 



1 

27" 



V^ (T) Z(r) = - Ric#(Z(r)) + -Z(t), (15) 



so that the last term in (1121) vanishes. □ 



Corollary 2 (Hessian of L; see [5j Corollary 7.40] for a similar statement). Let Z be a vector 
field along 7 of the form dHJ) and £ := Z(t{) © Z(^j) £ T^) (M x A/) Then 

Hess 9(ri)ffi9(T2) £| (a , |n . yi ^ (£,£) < - P V^H(j(r),Z(r))dT + ^f 



T=T\ 



2v^Ric 9(T) (Z(r),Z(r))|^. (16) 



Proof. Let V : (— e, e) x [n, T2] — )• M be any variation of 7 with variation field Z and such that 
V^( n )5(0, ri) and V^( T2 )5'(0,r2) vanish. Since 



ReSS 9(Tl)®9(T2) L \ {XiTV! y >n) (£>0 < ^2 



£(r s ), 

s=0 



the claim follows from Corollary [TJ D 

Let now Z* (i = 1, . . . , d) be space-time parallel fields along 7 satisfying Z* (n ) = uej 
(and consequently Z*(j2) = ve*), and Zj(r) := yr/tZ*{r) (so that £j = Zj(ri) © Zj(t2)). 
In order to estimate X^j=l ^- ess g(n)©3(T-2) -^L -1 (£«>£«) using Corollary [2] we will compute 

^ i=1 i?(7(r), Zj(r)) in the following (see [5j Section 7.5.3] for a similar argument). Set I\, I2 



and I3 by 



h:=-2Y,^{HT),Zi{T)% 



i=l 



a 

h ■= J2 [ ~ Hess i2(Zi(r), Z,(r)) + 2 |Ric # (Z 4 (r))|< 



i=l 



Ric(Z i (r),Z i (r))-2Rm(Z i (r),7(r),7(r),Z i (r)) 



7 3 := 4]T [(V Zi(r) Ric)(Z 4 (r),7(r)) - (V^ (t) Ric)(Z 4 (r), Z,(r)) 



i=l 



Then X^i=i ^(^(t), Zi{t)) = I\ + I2 + I3 holds. By a direct computation, 



I 2 = - t (-AR( 7 (t)) + 2 |Ric| 2 ( 7 (r)) - iu( 7 (r)) + 2Ric( 7 (r), 7 (r)) 
The contracted Bianchi identity divRic = ^Vi? [T3] Lemma 7.7] yields 

J 3 = y ((divRic)( 7 (r)) - (V^ (r)j R)( 7 (r))) = -^(V^ r) i?)( 7 (r)). 



For Ji, we have 

d 

/ 1 = - 2 £ 
i=\ 



(17) 



(18) 



dr 



(mc(Zi(T),Zi(T))) - {Vj {T) mc)(Z i (T),Z i (r))-2mc(Vj {T) Z i (r),Z i (r)) 



"2^ (^( 7 (r))) +2^ (T)J R( 7 (r)) + 4^Ric(V^ {T) Z J (T),^(r)) 

t=l 



dr Vt 



t 



2t / 1 
T Vr 



i?( 7 (r)) + ^( 7 (T))J+4^Ric(V i(r) Z i (r),^(r)). 



(19) 



Since Zj satisfies ([To]) . 
d 



4^Ric(V^ (T) Z i (r),Z i (T)) = 4^Ric(-Ric#(Z i (T)) + ^Z J (r),Z i (r)) 



i=l 



-^(2|Ric| 2 ( 7 (r))-ifi( 7 (r)) 



By substituting flU} into ([T9]) . 



7 i = -T (^(7(r))+2|Ric| 2 ( 7 (r)) 



Hence, by combining pi]). ([15]) and ([17]) . 



(20) 



(21) 



£ H(7(r), S*(t)) = I ( - 2^( 7 (r)) - 2 |Ric| 2 ( 7 (r)) - Ai2( 7 (r)) 
i=l 

-- J R( 7 (r)) + 2Ric(7(r),7(T))-2(V^ (T) ^)(7(r)) 
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Inserting this into (|16p we obtain 
d 

Yl HeSS <?(n)ffi9(r 2 ) L \( x , T1 ;y,r 2 ) (&> &) 



< ^ y ^ 3/2 ( v (t(t)) + 2 |Ric|2 (7(r)) + AjR ^( r )) 



T 



+ ^J2( 7 (t)) - 2 Ric(7(r), 7(r)) + 2(V^ t)j R)( 7 (t))) dr 



+ t 

d^ T - 



2r 3/2 

-^fl(7(T)) 



T=T1 



/ 



+ 7 / 2 r3/2 f 2 l Ric ! 2 (7(r)) + A J R( 7 (r)) - --R(7(r)) - 2Ric(7(r),7(r)))dr 
which completes the proof of Proposition [TJ 

5 Coupling via approximation by geodesic random walks 

To avoid a technical difficulty coming from singularity of L on the £-cut locus, we provide an 
alternative way to constructing a coupling of Brownian motions by space-time parallel transport. 
In this section, we first define a coupling of geometric random walks which approximate g(r)- 
Brownian motion. Next, in order to provide a local uniform control of error terms coming from 
our discretization, we study several estimates of geometric quantities in subsection I5.ll Those 
are obtained as a small modification of existing arguments in l23l [25] . The £-cut locus is also 
reviewed and studied there. Finally, we will establish an analogue of arguments in section [3] for 
coupled geodesic random walks to complete the proof of Theorem [2J 

Let us take a family of minimal /^-geodesies {'Jxy' 2 I ^i ^ r i < r 2 < ^2,x,y G M} so that a 
map (x, r\\ y, T2) *->• IxV 2 is measurable. The existence of such a family of minimal £-geodesics 
can be shown in a similar way as discussed in the proof of |15[ Proposition 2.6] since the family 
of minimal /^-geodesies with fixed endpoints is compact (cf. the proof of Lemma 7.27]). For 
each r G [ri,T], take a measurable section & T > of <7(r)-orthonormal frame bundle 9 ^ T '{M) of 
M. For x,y G M and t±,T2 G [tu^I with n < T2, let us define <&i(x, ri; y, T2) G F{M) for 
» = 1,2 by 

$i(x,ri;2/,r 2 ) := $ (ti) (x), 

*2(a?,ri;y,7i) := m^ T2 o $( ri )(x), 

where fnQ-S 2 is as given in Definition [2J Let us take a family of Revalued i.i.d. random variables 
(A n ) ng N which are uniformly distributed on a unit ball centered at origin. We denote the 
(Riemannian) exponential map with respect to g(r) at x G M by expi" . In what follows, 
we define a coupled geodesic random walk Xf = (X^ lt ,Y^ 2t ) with scale parameter e > and 
initial condition ~K e s = (x±,yi) inductively. First we set (X| lS ,Y? s ) := (xi,y%). For simplicity 
of notations, we set t n := (s + s 2 n) A (T/fj}). After we defined (X|) te r at i, we extend it to 
(Xf) te[s , Wl ] by 



C := Vd + 2$ l (X £ fltn ,f 1 t n ;Yl tn ,f 2 t n )\ n+1 , i = 1,2, 

ye ._ PYrl ( f i 4 «) / ^~^ A^FxC 1 ) 

y r 2 t — expy^ ^ £ V^-r 2 A n+1 j 
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We can (and we will) extend the definition of X%. for r G [Tfi/f2,T] in the same way. As in 
section [3] X^ t does not move when f\ = 0. Note that \/d + 2 is a normalization factor in the 
sense Cov(y/d + 2A n ) = Id. Let us equip path spaces C([a, b] — > M) or C([a, 6] — > M x M) with 
the uniform convergence topology induced from g(T). Here the interval [a, b] will be chosen 
appropriately in each context. By (|23|) which we will review below, different choices of a metric 
<?(r) from g(T) always induce the same topology on path spaces. As shown in jTTj , (X!f) Te \ flSt T} 
and (y r £ ) TG [ f2S pn converge in law to <7(r)-Brownian motions (^t) t g[tis,T] an d 0^r) T &[f 2 s,T] on M 
with initial conditions Xf lS = x±, Y f2S = y\ respectively as e — > (when f\ > 0). As a result, 
X £ is tight and hence there is a convergent subsequence of X e . We fix such a subsequence 
and use the same symbol (X e ) e for simplicity of notations. We denote the limit in law of X e 
as e — > by X t = (Xf lt ,Y f2 t). Recall that, in this paper, g(r)-Brownian motion means a 
time-inhomogeneous diffusion process associated with A 5 ( T ) instead of A g ( T )/2. 

Remark 8. We explain the reason why our alternative construction works efficiently to avoid 
the problem arising from singularity of L. To make it clear, we begin with observing the essence 
of difficulties in the SDE approach we used in section [3j Recall that our argument is based on the 
Ito formula. Hence non-differentiability of L at the £-cut locus causes the technical difficulty. 
One possible strategy is to extend the Ito formula for ^-distance. Since £-cut locus is sufficiently 
thin, we can expect that the totality of times when our coupled particles stay there has measure 
zero. In addition, as that of Riemannian cut locus, the presence of £-cut locus would work to 
decrease the ^-distance between coupled particles. Thus one might think it possible to extend 
Ito formula for /^-distance to the one involving a "local time at the C-cut locus" . If we succeed 
in doing so, we will obtain a differential inequality which implies the supermartingale property 
by neglecting this additional term since it should be nonpositive. 

Instead of completing the above strategy, our alternative approach in this section directly 
provides a difference inequality without extracting the additional "local time" term. By dividing 
a minimal ^-geodesic into two pieces, we can obtain a "difference inequality" of ^-distance even 
when the pair of endpoints belongs to the £-cut locus (see Lemma [3|). In order to employ such 
an inequality, it is more suitable to work with discrete time processes. 

5.1 Preliminaries on the geometry of ^-functional 

Recall that we assume that M has bounded curvature, so that there is a constant Co < oo such 
that 

, , mp r ^ I Rm Iff(r) ( x ) v I Ric l ff (r) 0*0 < Co- (22) 

On the basis of (j22[) . we have a comparison of Riemannian metrics at different times. That is, 
for n < r 2 , 

e _ 2Co ( r2 _ ri ) 5(r2) < g{ji) < e 2C (T 2 -T l)g{T2 y (23) 

Let p g ( T ) be the distance function on M at time r. Note that a similar comparison between 
p g ( T1 ) and Pg(j 2 ) follows from (f23|) . We also obtain the following bounds for L from ([22"]) and 
([23]) . Let 7 : [ti,T2] — > M be a minimal £-geodesic. Then, for r € [tl,T2], 

e -2C (r 2 -Ti) 2 3/ 2 3 / 2 

--Pg(T)(l(ri),7(r)) --dC (T 2 ' - t^ ) < L(j(t 1 ),t 1 ;j(t 2 ),t 2 ) 



2Jn-Jn ry ^' XIK ^" v " 3 



2Co(r 2 -ri) 9 o/o o/o 

-P 9 (T)(7(n),7(T 2 )) 2 + -dC (r 3 2 /2 - rl 12 ) (24) 



(see [5] Lemma 7.13] and [231 Proposition B.2]). The same estimate holds for p g (T){l{ T )-,l{ T 2)) 2 
instead of Pg<T){l{ T i)il{ T )) 2 ■ Taking the fact that /^-functional is not invariant under re- 
parametrization of curves into account, we will introduce an estimate for the velocity of the 
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minimal /^-geodesic 7. By a similar argument as in [5] Lemma 7.13 (ii)], there exists r* £ [ti,t 2 ] 
such that 

r*|7(OG(r.) < 2(^ 1 -^F 1 ) (^(7(Ti),r 1 ; 7 (T 2 ),r 2 ) + ^(rf 2 -r 1 3/2 )) . (25) 

Suppose r 2 < T. Then, as shown in [21] (see [6] also), ([22]) yields that there is a constant 
C(d) > depending only on d such that 



C(rf)Co 

t<t 2 ,x£M (T — r 2 ) A Cq 



sup |V Rm| 9(T) (*) < - ^ ^°_ =: cj. (26 ) 



By virtue of ([26]) . there exists a constant c\, C\ > which depends on Co, C and T such that 
for all t[, t' 2 £ [ti, r 2 ] with r{ < t' 2 , 



t 2 \i(4)\; {ri) < PiTi|7(Ti)|; w) + C U (27) 

r{ W^l^^c^ 17(^)1^,) +Ci. (28) 



The first inequality in (|27|) can be shown similarly as Lemma 7.24]. It is due to a differential 
inequality based on the £-geodesic equation (J7|) which provides an upper bound of 5 r (r|7(r)P/ r J. 
By considering a lower bound of the same quantity instead, we obtain the second inequality (j28|) 
in a similar way. Combining (|27p and (|28p with (|25p and (|24p . we can take constants c 2 > and 
C 2 > depending on Co, c\, C±, n and r 2 such that 

r|7(r)| 2 (r) < c 2 p 9(T) ( 7 (r 1 ),7(r 2 )) 2 + C 2 (29) 

for ri < r < r 2 . Though c 2 and C 2 depends on t\ and r 2 , it is easy to see that c 2 and C 2 are 
uniformly bounded above as long as r 2 — t\ and T — r 2 is uniformly away from 0. 

Let us recall the definition and some properties of £-cut locus according to [5], I23|, [25] . Given 
r, t' £ [0, T) with t < t' and x £ M, we define /^-exponential map £ T ,r' ex Px : T X M — > M 
by C T y exp x (Z) = 7(1"'), where 7 is a unique £-geodesic from (r, x) with the initial condition 
lim T /| r Vr7(r') = Z. Note that we can extend the domain of C- geodesic 7 to the interval [r, T) 
by using ([27]) (see j5[ Lemma 7.25]). Set 



n(x,T 1 ;T 2 ) :=<Z£T X M 



7 : [n,T 2 ] ->• M given by 7(7-) := £ n]T exp x (Z) 
is a minimal £-geodesic 

r(x,r; Z) := sup{r G (n,T) |ZG $7(x,ti;t)} . 



Based on these notations, we define the £-cut locus £Cut by 
£Cut := { (x,T 1 ;y,T 2 ) 



xeM, n e [fi,T), 

y = £>ti,t 2 e W x (Z) for some Z S T X M, 

r 2 = f(x,n;Z) e (n,r) 



As remarked in [5] l23l l2"5], £Cut is a union of two different kinds of sets. The first one consists of 
(x, n; y, r 2 ) such that there exists more than one minimal £-geodesics joining (x, n) and (y, r 2 ). 
The second consists of (x,Ti;y,T 2 ) such that (y, r 2 ) is conjugate to (x,ti) along a minimal £- 
geodesic with respect to £-Jacobi field. Note that we can define exponential map in the reverse 
direction in r. By using this reverse exponential map, "reversed £-cut locus" is defined and it 
is identified with £Cut by virtue of the above characterization of £Cut. 
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5.2 Proof of Theorem H 



For the proof of our main theorem based on discrete approximation, we will follow a similar way 
as in previous studies in this direction (see [101 E] and references therein) . Our first task is to 
show a difference inequality of A(£, Xf ) in LemmaEJ We begin with introducing some notations. 
Set 7„ := 7xe"' T2 ™ an d let us define a vector field A„ +1 along j n by A n +i(r) = ^/r/t^A* , 1 (r), 

where A* +1 is a space-time parallel vector field along 7 n with initial condition A* +1 (fi£ n ) = X^+v 
Let us define random variables Q n and S n as follows: 



;t 



Cn+i := V2r (A^ +1 (r),7 n (r)) s(r) 
1 



r=fit n 



Sn+i := r r 3 / 2 (i? g(T )(7n(r)) - |7n(r)| 2 (r) ) 



+ 



T-2,t n 
T =f 1 t n 

t 



2/FRic 9(r) (AT +1 (r), Aj, +1 (r)) 



T-iX n 
T=f\t r , 



^H(j(T),Xl +l (r))dr). 



Tltn 
Tlt„ 

Here if is as given in ([13j). For M C M, we define ctm : C([s, T/f 2 ] -> M x M) -> [0, oo) by 

a Mo (w,w) : = inf {i > s | lo^t ^ M or tu^ t ^ M } . 

For simplicity of notations, <tm (X £ ) and <tm (X) are denoted by a E M and <XjL respectively. As 
shown in [TT], for any r\ > 0, we can take a compact set Mq C M such that lim^o P[o"m <T]<r) 
holds (cf. [12]). 

Lemma 3. -Lei Mo CMie compact. Then there exist a family of random variables (Qn)neN,e>o 
and a family of deterministic constants (8(e)) e> o with lim e _j,o <5( £ ) = satisfying 



such that 



n; tn«^ o A(T/f 2 ) 

A(t n +i,Xf n+i ) < A(i n ,Xf n ) +eCn+i + £ Sn+i + Q^+i- 



(30) 



(31) 



Proof When (X £ (fit n ), T\t n ; Y £ (f2t n ) , f2t n ) ^ £Cut, the inequality (pJT]) follows from the Taylor 
expansion with the error term Q £ n+ \ = o{e 2 ). Indeed, the first variation formula ([5j Lemma 7.15] 
cf. dHJ) produces eCn+i and Corollary [2] together with @ implies the bound e 2 S n+ i of the second 
order term. To include the case (X £ (fit n ),fit n ;Y e (f2t n ),f2t n ) £ £Cut as well as to obtain a 
uniform bound ()30p . we extend this argument. Set t* := {f\ + f2)t n j1. Then we can show 

( X fl tn > n f " J 7n (t* ) , T * ) £ £Cut , 

(7 n (r*),^;X| 2tn ,f 2 t n )^£Cut 

since minimal £-geodesics with these pair of endpoints can be extended with keeping its mini- 
mality (cf. see [3 Section 7.8] and [25]). Set x* n+l = exp (r y M ( y/fi + f 2 a! . i (r*) ) . The triangle 
inequality for L yields 

A ( i n,X^ n ) = L(X| lin ,fit n ;7 n «),r*) + L(-y n (T*),T*;X e . 2tn ,f 2 t n ), 
A(i n+ i,Xf J < L (X| lin+i ,fit n+ i;< +1 ,r* +1 ) + L (< +1 , r* +1 ; X^ tn+1 ,f 2 t n+1 ) . 
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Hence 



A(t n+ i,Xf n+i ) - A(t n ,,XfJ < [L [X e fitn+v Tit n+ i;x* n+1 ,T^ +1 J - L(X| lfn ,rit n ;7„(r*),r, 

and the desired inequality with Q e n = o(e 2 ) holds by applying the Taylor expansion to each term 
on the right hand side of the above inequality. 

We turn to showing the claimed control (|3Up of the error term Q e n . Take M\ D Mq compact 
such that every minimal £-geodesic joining (x, fit) and (y, f 2 £) is included in Mi if x, y £ Md 
and t € [s, T/f 2 ] . Indeed, such Mi exists since we have the lower bound of L in (|24|) and L is 
continuous. Let us define a set A by 



A:~- 



((ri,x),(T3,z),(r 2 ,y)) G ([fi,T] x Mi 



x,y € M , 

72 - Ti > (f 2 - fi)s, 

T 3 = (n + T 2 )/2, 
L(x,ti;z,t 3 ) + L(z,T 3 ;y,T 2 ) 
= L(x,Ti;y,T 2 ) 



> . 



Note that A is compact. Let vti,7T2 : A — >• ([ri,T] x M±) 2 be defined by 

7Ti((ti,x),(t3,z),(t2,?/)) := ((ri,x), (73,2:)), 
7r 2 ((ri, x), (r 3 , z), (r 2 , y)) := ((r 3 , 0), (r 2 , J/)). 

Then 7Ti(>l) and tt 2 (>1) are compact and iri(A) n £Cut = for i = 1,2. The second asser- 
tion comes from the fact that (2,73) is on a minimal £-geodesic joining (x,ri) and (y,r 2 ) for 
((x,ri), (z, T3), (y,r 2 )) G A. Note that £Cut is closed (see [23]; though they assumed M to be 
compact, an extension to the non-compact case is straightforward). Thus we can take relatively 
compact open sets Gi,G 2 C [fi,T] X M such that TTi(H) C Gj and Gi n £Cut = for i = 1,2. 
Then the Taylor expansion we discussed above can be done on G\ or G 2 for sufficiently small e. 
Recall that L is smooth outside of £Cut (see [5]). Thus the convergence e~ 2 Q n (e) — > as e — > 
is uniform in n and independent of X| as long as t n < a e M A (r/f 2 ). Since the cardinality of 
{n J t n < a £ M A (T/f 2 )} is of order at most e~ 2 , the assertion ([30]) holds. □ 

We next establish the corresponding difference inequality for ©(£, Xf) (Corollary [3|). For 
that, we show the following auxiliary lemma. 

Lemma 4. Let T < T. 

(i) There exist deterministic constants c 2 > and C 2 > such that 

|Cn|<cV ff (T)(Xt_ 1 ) + C 2 , |A(t n ,X|J| <c 2/VT) (X?J 2 + C 2 

('mJ TaA;e Mo C M compact. Then there is a constant R = R(Tq,Mq) > such that |S n | < i? 
ZioWs i/t n < o-| /o A (T /f 2 ). 

Proof. By the definition of £ n , we have 



|Cn| < V 2 ( d + 2)*n-l (n|7n-l(ri* n -l)| a (T l t n _ 1 ) +75|7n-l(75*n-l)| ff (fl,t,»_i)) ■ 



Thus the desired bound for |^ n | follows from (|29p and ()23p . Similarly, the estimate for A(i n , Xf ) 
follows from (|24[) and (|23|) . For the assertion (ii), we deal with the integral involving H in the 
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definition of £ n . Note that every tensor field appeared in the definition of H is continuous. As 
in the proof of Lemma [31 take Mi D Mq compact such that every minimal £-geodesic joining 
(x, fit) and (y, f 2 t) is included in Mi if x, y G Mo and £ G [s, T/f^]. Since Xf n _ 1 G Mo x Mo holds 
on the event {£„ < crf /o A(To/f2)}, the upper bound (|29|) of V^ItC 7 ")! implies that H (j n (r) , Z (t)) 
is uniformly bounded for any vector field ■Z'(t) along 7 n of the form Z(t) = \Jt /t n Z*(r) with a 
space-time parallel vector field Z*(t) satisfying \Z*{r)\ g ^ < 1. 

This fact yields an expected bound for the integral. For any other terms in the definition of 
E n , we can estimate them as in the assertion (i). □ 

By virtue of Lemma 01 A(i n ,X| ), C, n and E n are uniformly bounded on the event {t n < 
a M ^ (^0/^2)} for To < T. Thus Lemma [3] yields the following: 

Corollary 3. Let T$ <T and M$ C M be a compact set. Then there exist a family of random 
variables (On)neN,e>o an d a family of deterministic constants (6(e)) £> o with lim E _>.o <f(e) = 
satisfying 

E &n< He) 

n; t n <(X e Uo A(To/f2) 

such that 



e(t n+1 ,xf n+1 )<e(t n ,xfj + 



e 



'■11 



ffi-y/n) A(t n , Xf J - 2e 2 d ( v 7 ^ - V^) 



+ 2e ( y/T2t n +i - ^Tit n+ ij Cn+1 + 2e ( y / T 2 t n+ i - yViVf-lJ S n+i 

(32) 



+ QL 



The term S n corresponds to the one dominating the bounded variation part of dA(t, Xt,Yf) in 
section[3l However, as a result of our discretization, we are no longer able to apply Proposition Q] 
directly to estimate S n itself. In this case, we can do it to the conditional expectation of E n 
instead. Set Q n := cr(Ai, . . . , A n ) and E n +i := E[S n+ i | Q n ]. Then, since each <!>,; is isometry and 
(d + 2)E[(A n ,ej)(A n , ej)] = Sij, Proposition Q] yields 



d 



1 






(33) 



In order to replace S n with S n in (|32p . we show the following: 

Lemma 5. Let to < Tq < T and Mq C M compact. For t G [^i,T], set Nf := sup{n G 
N I f 2 (s + e 2 n) < t}. Then, for r] > 0, 



lim sup I 

£-»0 



sup 

Nf<N<N£ 

H) - - ~0 
. tN<<r e M 



N 



/ u ytnV-'n ~ E r 



71=1 



> £ 77 



0. 



Proof. Lemma 2] ensures that S n and £ n is bounded as long as n < N™ and t n < a e M . Thus 
X^n=i V^n(^n ~~ ^n) is a ^Ar-local martingale. Hence the Doob inequality implies 



lim sup ] 

e->0 



sup 

0<N<NI 
tN<T e Mo 



N 

E 

n=l 



't n (E r 



> e 77 



(34) 
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for t £ [ti,To]. Since we have 

N 

/ j \~b n (lJ n — 2_l n ) 



sup 
<j 

r M 



Nf <N<NZ 



?i=JVf 



> e i] 



C ^ sup 

0<N<NL 



t N <a\ 



<<> 



A 1 



/ vin(S n — T, n/ 



ra=l 



> 



£ 2 7] 



> u < 



sup 

o<N<m 



N 



/ J yt n (lj n 2->n, 



n=l 



> 



£ 2 7] 



tjV<cr 



M (j 



the assertion follows from ([34 

As a final preparation, we show the following auxiliary lemma. 
Lemma 6. There exists C3 > such that 



□ 



E 



sup |e(t,x t )| 2 

s<t<T/f 2 



<c 3 . 



Proof. By virtue of ()24p . 0(i, Xt) is bounded from below uniformly in t £ [s, T/t2[. In addition, 
there is a constant c, C > such that 

G(t,X t ) < cp 9(T) (X t ) 2 + C 

holds for t £ [s,T/f2J. Take a reference point o £ M. Then we have 

/5 9 (T)(Xt) < e CoT (p g{Tlt) (o,X Tlt ) + p g{T2t) (o,Y T2t )) . 

Thus the proof can be reduced to the following claim: 



E 



sup p g (T 2 t)(o,Y f2 t) 

f 2 s<t<T 



< 00. 



(35) 



Indeed, a similar bound for X flt follows in the same way. As shown in [12], (p g ^(o, Yt))[f 2 s,T] 
is dominated from above by a Bessel process (of dimension d) plus a constant. Thus ([35]) 
easily follows from the Burkholder inequality for the fourth moment of a Euclidean Brownian 
motion. □ 

Proof of Theorem® First we remark that the map (x,y) *— > (X a ,Y a ) is obviously measurable. 
Thus, we obtain the same measurability for (X, Y). The integrability of 0(£,Xj) follows from 
Lemma El We will show the supermartingale property in the sequel. For s < s± < ■ ■ ■ < s m < 
t' < t < T and /1, . . . , f m £ C C (M x M -)■ R) with < fj < 1, Set F(w) := UT=i /iK) 
for w £ C([s,T/f2] — > M x M). Take r] > arbitrarily and choose a relatively compact open 
set M C M so that P[o^ /o <*]<?? holds. Note that limsup £ ^ P[cr| fo < t] < r] also holds 
since {w \ o~m ( w ) ^ t} is closed. It suffices to show that there is a constant C > which is 
independent of 77 and Mq such that, 



E 



'(@(tAa° Mo ,X tA<T o Mo ) - @(t'Aa Mo ,X t , Aff o Mo )) FfX.^)] < C0? 

holds. In fact, once we have shown (|36p . then Lemma [6] yields 

E[(G(t,X t )-G( S ,X s ))F(X)]<0 



(36) 
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since a° M —> oo almost surely as Mo f M. 

Take / G C C (M x M) such that < / < 1 and f\ v = 1, where [/ C M x M is a open set 
containing Mo x ^o- Then, by virtue of Lemma [6] and the choice of Mq, 

E[(@(tA^ /o ,X tAff o /o )-e(t'AaO /o ,X t , ACT o /o ))F(X. A(T o /o 

< e [(e(t,Xt) - e(i',x t ,)) /(x t )/(x t o J F(x) ; ^ /o >t]+ i&J 2 ^. (37) 

For it G [s, T/f2], let us define \u\ e by 

|_uj £ := sup{s + e 2 n | n G NU{0},1 + e 2 n < u}. 
Then, since {w \ (Jm {w) > t} is open, 



e [(e(t,Xt) - e(t',x t o) /(x t )/(Xt/)F(x) ; <, > t] 

< liminf E [(G(t,Xf) - 6(i',Xf,)) /(X?)/(X£)F(X e ) ; a^ > t] 



lim inf E 
e-s-0 



6(LtJ £ ,Xf tJe ) - e(Lt'J e ,Xf t , Je ) /(Xf tJ J/(X e Lt , Je )F(X £ ) ; a% Io > t . (38) 



Here the last inequality follows from the continuity of and /. Set a £ M := \& e M \ £ + e 2 . Note 
that {a e Mo = t n } G Q n for all n G N U {0}. Then 



E 



9(LtJ„Xf ij£ ) - 8(|/J e ,Xfcjj) /(X e LtJ J/(X^ Je )F(X £ ) ; a% Io > t 



< E 



Q([t\ e Aa% Io ,Xl tiEAah )-G(Lt / J,A^ /0 ,X| i , JeA - i/ ) F(X 



■A^ 



+ 2E 



sup |e(«,x£)/(x£ 

s<u<T/t 2 



em2 



1/2 



[a E Mo < t] 1/2 - (39) 



Since a function w i-> sup 1<u<T/f2 \@(u,w u )f(vr u )\ on C([s,T/t2] — )■ M x M) is bounded and 
continuous, we have 



lim sup E 

£^0 



sup |0(«,X£)/(X£) 

s<n<T/-f 2 



1/2 



[^ /0 < t) 1 ' 2 < cfV* 



(40) 



By combining ([38]) , ([39]) and ([40]) with ([37]) , the proof of ([36]) is reduced to show the following 
estimate: 



lim sup E 

£^0 



(9(L*J e A4,X| tJeA ^) - @([t'\ £ Aa^.X^^)) F(X £ A - |/() )] < C^7- (41) 



Take N\ and A^ so that t^ = \t'\ e Aaj^ and £at 2 = [t\ £ Aa e M hold. Let .E^ be an event defined 
by 

N 2 
ra=JVi 



-En 



< 



2e 2 ( v / ^-\/n) 



On £^, an iteration of ([32]) together with ([33]) yields 



N-2 



@(t N2 ,X tN2 ) - @(t Nl ,X tNl ) < e ^ (V^n ~ Vntn) Cn + 2y/V 

n=Ni+l 
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for sufficiently small e. In addition, Lemma [5] yields limsup^^Q P[E£] = 0. By applying LemmaE] 
with Tq = fit to an iteration of (|32p . we obtain a constant C > satisfying 



<C 



N 2 
n=Ni+l 



uniformly in sufficiently small e > 0. Since F(X. £ A ~ ) is Qn 1 -measurable, we obtain 



E 



(G^X? )-e(^,X? ))F(X £ A -, ) <C¥[E^] i/z + 2^}. 



1/2 



Hence (|4ip holds with C = 1/2 and the proof is completed. □ 
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